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We discuss the consequences of the accretion of dark matter (DM) particles on compact stars 
such as white dwarfs and neutron stars. We show that in large regions of the DM parameter space, 
these objects are sensitive probes of the presence of DM, and can be used to set constraints both 
on the DM density and on the physical properties of DM particles. 



I. INTRODUCTION 

Observations continue to provide overwhelming ev- 
idence for a dark component of matter in the Uni- 
verse [H [2} . There are a number of different candidates 
for this dark matter (DM), for example, weakly interact- 
ing massive particles (WIMPs) are attractive candidates 
because their weak scale annihilation cross-section natu- 
rally gives rise to densities in today's universe comparable 
with what is needed cosmologically. They are also inter- 
esting candidates in as much as they may be probed by 
the increasingly accurate direct detection experiments in 
development such as the current leader XENON10 [|[. 
Furthermore, it is possible that WIMPs could be pro- 
duced at the soon to be operating Large Hadron Collider 
(LHC) at CERN. 

There also exist other candidates which are more 
weakly coupled to standard model particles. The axion is 
an example of a very feebly coupled cold DM candidate, 
whereas the gravitino or the sterile neutrino can be can- 
didates for either cold or warm DM, depending on their 
production mechanisms and masses and are equally dif- 
ficult to detect (see Refs. [H, [1] and references therein). 
Another set of DM particles which have been suggested 
are superheavy particles known as WIMPzillas which are 
produced non-thermally at the end of inflation 0, [f| . 

It has long been realised that a finite DM-nucleon cross 
section would result in scattering between DM and the 
gas in stars, knocking the DM particles into orbits which 
result in their subsequent capture [1, 0, H, Q • Such cap- 
tured particles would "sink" to the centre of stars where 
they would annihilate with themselves. Indeed neutrinos 
from such annihilations at the centre of the sun will be 
searched for at neutrino telescopes such as ICECUBE, 
currently under construction in Antarctica [lj| HH ■ 

The amount of DM in the solar system is such that 
while one can expect these products of WIMP annihila- 
tion in the centre of the Sun, the total injected energy 
rate associated with this annihilation is a tiny fraction of 
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the luminosity in the sun due to nuclear burning. We do 
not therefore expect that the presence of these WIMPs 
will affect the way that the Sun itself burns, although 
future helioseismic observations may probe an interest- 
ing region of the WIMP parameter space[l2|. Here, we 
focus on situations where the capture of WIMPs may sig- 
nificantly affect the strcture and/or the evolution of the 
astronomical objects onto which they are captured. 

N-body simulations of DM halos suggest that the DM 
density close to the galactic centre should be larger by 
several orders of magnitude than the density in the so- 
lar system (see e.g. the recent "Via Lactea" simulation 
of Diemand et al. [IH and references therein) . Stars in 
those regions therefore may experience the accretion of 
a large amount of DM which could change the internal 
energy exchange mechanisms, e.g. changing the temper- 
ature gradient in convective stellar cores [141 ]. Further- 
more, if the density of DM is extremely high, the energy 
injection rate due to DM annihilations could compete 
with the usual nuclear energy production rate (l5l lid ] . 
As for the Earth, the requirement that the injected en- 
ergy does not exceed the measured Earth's heat flow, has 
been used to set interesting constraints on the scattering 
cross-section of DM particles off nuclei [TtJ • 

The capture of DM onto stars is proportional to the 
number of nucleons in the star times the escape veloc- 
ity, therefore compact objects such as white dwarfs or 
neutron stars are ideal targets for searches aimed at de- 
tecting the effects of DM accretion. Both classes of de- 
generate compact objects are unfortunately also usually 
rather hot, so that any increase in temperature due to the 
accretion of WIMPs would be difficult to detect. In the 
past decade, however, great advances have been made in 
the detection of cool white dwarfs in globular clusters, 
where very old stars with temperatures of less than 4000 
K have been observed [HI]. It would seem therefore to 
be a good place to look for the effects of WIMP accre- 
tion and indeed we will show that the constraints which 
can be obtained are very interesting, despite being model 
dependent. 

As for neutron stars, we will see that it is rather diffi- 
cult to detect the heating produced by the annihilation 
of WIMPs captured in their cores, since the total cross 
section for capture is limited by the relatively small size 
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of the star. The accretion of DM particles which are not 
their own anti-particles or do not annihilate with them- 
selves very quickly is however constrained, since in this 
case DM particles (if they are more massive than nu- 
cleons and have couplings to them) would accumlate at 
the centre of the neutron star, eventually becoming the 
dominant source of gravitational potential in that central 
region. As we shall see, if the density of DM in the core 
were to increase more rapidly than its self annihilation 
rate, then the accumulation could lead to a gravitational 
instability and the collapse of the DM cloud. 

The paper is organised as follows: in the next section 
we will discuss the accretion of DM particles onto com- 
pact objects in generality. In section IIIII we will then 
develop a model for the baryons and the DM in globular 
clusters, specifically M4, where observations of cool stars 
have been made close to the centre of the cluster where 
the DM, if there at all, is expected to be densest. We 
argue that if there is DM in globular clusters, then its 
density at the radius we are interested in is an extremely 
weak function of the total mass within the cluster. We 
work out the accretion rate of DM onto white dwarfs and 
estimate how it would change their luminosity, finally 
comparing with the actual data to place constraints on 
the combination of DM density and WIMP-nucleon cross 
section. 

Then in section ITVl we will describe the various possi- 
ble effects upon neutron stars due to the accumulation 
of DM. We calculate the expected luminosity of galac- 
tic neutron stars at different distances from the centre 
of the milky way, before considering the interesting ef- 
fects which could happen inside neutron stars due to the 
accretion of a DM species with a small self annihilation 
cross section. 



II. ACCRETION OF DM ONTO GRAVITATING 
OBJECTS 

The accretion of WIMPs onto stars has been stud- 
ied for many years now [(J 0, H, H|. The capture rate 
Gamma c which we adopt interpolates between the two 
different capture rates defined for optically thin and thick 
bodies in Ref. 19] 
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where pdm is the density of DM round the star, md m is 
the mass of the DM particles and v is the average DM 
velocity. v esc is the escape velocity of the astronomical 
body in question. 

The effective cross section cr e ff is given by by one of 
two things, normally it is the sum of the cross sections 
upon the individual nuclei in the star, including the co- 
herence factor which means that the spin independent 
cross section of a nucleus of element i with A^ nucleons 
will be Af times the cross section a s t of an individual 



proton or neutron (see e.g. [20|, l2l|). However, this to- 
tal cross section cannot be larger than the geometrical 
size of the star, so the effective cross section is given by 
whichever is smallest 
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where xi is the mass fraction of element i, M* and i?* 
are the mass and radii of the star and m p is the proton 
mass. For the capture of WIMPs onto neutron stars, 
one expects that the geometrical cross section may be 
less than the sum of the individual cross sections. The 
cross section for scattering between the DM particle and 
nucleons would have to be less than 10~ 45 cm 2 in order 
for the sum of the individual cross sections in the star to 
be less than the geometrical cross section of the star. 

Once DM is captured by a star or a compact object, it 
will form a thermal distribution inside the star of char- 
acteristic radius 191 
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where p c and T c is the central density and temperature 
and mdm is the mass of the DM particle. Typical values 
of this thermal radius for WIMP DM will be around a 
metre for a neutron star or a kilometre for a cool white 
dwarf. Having been concentrated in the centre of the 
star, the DM will annihilate with itself, the equation for 
the evolution of number density in the star over time is 
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where the annihilation rate T a is given by 
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which contains the DM self annihilation cross section 
(a ann v) and the total number of DM particles in the star 
N. When r c = T a the capture and annihilation rates will 
be equal and any DM which accretes onto the star will 
be instantly converted into additional luminosity. The 
timescale for this steady state to be reached is given by 
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In the next section we will go on to calculate the effect 
of DM accretion onto cool white dwarfs. 



III. WHITE DWARFS IN GLOBULAR 
CLUSTERS 

The presence of DM in globular clusters is a contro- 
versial issue, with early work by Peebles [23] suggesting 
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that globular clusters form along with dark halos. Since 
then there are those that argue that even if the clus- 
ters did start with a halo of DM that it will have been 
tidally disrupted since then due to interactions with the 
host galaxy [24(. However, other studies suggest that 
the cores of DM halos in globular clusters will survive 
successive tidal interactions with the host galaxy [25|. 
The combination of weakly interacting DM and old cool 
white dwarfs therefore seems to present a real possibility 
of testing the properties of DM. 

The distribution of stars in clusters can be fit well by 
solutions of the collisionless Liouville equation for a given 
velocity distribution [26[. The cluster profiles are then 
determined by the core radius r c and the tidal radius rt 
plus an overall normalisation. The relevant parameters 
for the globular cluster M4 are a core radius of r c = 
0.83' in arc-minutes and a concentration parameter c = 
log(rt/r c ) = 1.59 [22]. This, combined with the distance 
to the cluster of 1.73 kpc [18|], gives r c — 0.417 pc. We 
use these parameters to model the density of stars in 
the globular cluster using the King model. The overall 
normalisation of density is obtained by integrating the 
King profile and setting the total mass to 1O 5 M , this 
figure obtained from the apparent magnitude and the 
distance which give a total luminosity of the cluster of 
L ~ 1.5 x 10 5 L Q . 

In order to estimate the size and the extent of the DM 
halo surrounding the globular cluster, we follow the pro- 
cedure of [25| . The average cosmological ratio between 
baryons and DM is ^lb/^DM ~ 0.2 [28l | while the fraction 
of baryons in the form of globular clusters is F = 0.0025 
[29l | . Mashchenko and Sills have therefore estimated that, 
having taken into account mass loss during stellar evo- 
lution (and other factors, see Ref. [2^| for details), the 
ratio of baryonic to non-baryonic mass is F — 0.0088. 
Here, we assume that the mass of the DM in the Glob- 
ular cluster is Mum = 10 7 M Q , and we will show below 
that our results are rather insensitive to the exact value 
of the baryonic fraction. 

We define the virial radius R V i r as being the radius 
within which the average density is bigger than the av- 
erage cosmological density by a factor A, called virial 
over 'density , for which we have adopted here the fitting 
form of Bryan and Norman [30T| 



A = 18tt 2 + 82(ft m - 1) - 39(O m - l) 2 



101 



leading to the relationship 
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where we have assumed a Hubble constant of 72 
s _1 Mpc 
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and Q m = 0.3. A virial radius at z = 
of more than 5 kpc does not really make any sense be- 
cause globular clusters are well embedded within the 
much larger halo of the milky way, but evidence from 
N-body simulations suggests that the inner part of DM 
halos survives tidal stripping [25| . and we will only use 
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FIG. 1: Density of baryons and DM in the central region of 
the globular cluster M4. The solid line shows the density of 
baryons as a function of radius. The dotted and dashed lines 
represent the density of DM for the normal NFW profile and 
for an adiabatically contracted profile, respectively. 



R V i r to derive the properties of the initial central DM 
distribution. Following [32|, [33j], we will assume a con- 
centration parameter for a halo of this mass to be c = 32 
(slightly different models give rise to different c values, 
but such changes will not change the results significantly) 
and we will assume a Navarro, Frenk and White profile 
of the form SI 
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with a = 171 pc and 

Mum 



Pc 



4ira 3 



ln(l + c) 



l + c 



(9) 



0.63M Q pc- A (10) 



which finishes our definition of the DM profile. The dif- 
ferent density components are shown in Fig. [TJ 

We will be interested in the density of DM close to 
the core radius r c of the globular cluster. As we claimed 
earlier, this is a very weak function of the total mass of 
the halo, Mdm- In fact p(r c ) changes by around a factor 
3 for different halo masses between 10 6 Mq and 10 8 Af Q . 
The reason for this weak dependence on the total mass is 
that the density of DM at some small fixed radius r -C a 
is proportional to p c a which in turn is proportional to 

m dm c2 /(M± + c) - c/(l + c)) and as M DM goes down, 
the halo profile concentration c increases, see for exam- 
ple [33j]. This will make our results less sensitive to the 
overall mass of DM. 
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The DM halo will be distorted by the presence of 
baryons, leading to an enhancement of the halo profile 
in the inner regions of the globular cluster [35|, Ha, H3| . 
At the same time, the heating of DM particles due to 
interactions with stars will tend to wipe out the effect of 
adiabatic contraction, over a timescale [3811 
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where e in the kinetic energy of a typical DM particle, p* 
and m* are the density and mass of stars, G is Newton's 
constant, v rrns is the rms velocity of the DM particles 
which we assume is equal to that of the stars and A ~ 
0.47V where N is the total number of stars below the 
radius in question. At the core radius of the cluster, the 
timescale for heating of DM is much less than the age 
of the universe, so one cannot expect to rely upon the 
density profiles plotted in figure Q] since globular cluster 
MA is nearly as old as the Universe. We will therefore 
assume a maximum density of 50 Mq pc -3 , which is the 
approximate density at the radius where the timescale 
for heating of DM becomes larger than the age of the 
Universe. 

A 0.5Mq white dwarf has a radius of approximately 
10 9 cm and will be composed of carbon or oxygen (the 
densities and observed temperatures of newly formed 
white dwarfs are inconsistent with the presence of hydro- 
gen or helium, see e.g. 22|). Since we wish to make as 



many conservative assumptions as possible in the calcula- 
tion, and because the white dwarfs we are looking at are 
formed in low metallicity environments, we will assume 
that the white dwarfs in question are composed entirely 
from carbon. We can then work out the spin indepen- 
dent cross section which would be required in order for 
the geometrical cross section to dominate over the sum of 
the individual nuclei cross section, which turns out to be 
around 10~ 41 cm 2 . Because the constraints on the spin 
independent WIMP-nucleon cross section from XENON 
are at the level of 10 -43 cm 2 , we will assume that a e ff is 
given by the sum of the individual nucleon cross sections 
rather than the geometrical cross section of the star. 

Other ingredients that we need for our calculations are 
the escape velocity of the white dwarf, which is trivially 
obtained, and also the velocity of the DM particles at 
the radius corresponding to r c . Studies of the velocities 
of stars in globular cluster M4 show typically a velocity 
dispersion of a few kilometres per second [39(. Equa- 
tion fl} shows us that the capture rate increases for low 
velocity DM, so to be conservative we choose the max- 
imum velocity possible at r c which corresponds to the 
escape velocity obtained by integrating the derivative of 
the potential of the King profile of stars (as in the Milky 
Way, the inner potential of M4 is completely dominated 
by baryons). This gives us a velocity for the DM of 20 
km s _1 , assuming that there is no black hole at the cen- 
tre of the cluster. The presence of such an object would 
anyway seem to be at odds with the observed velocity 
dispersion of stars. 



The timescale for the equilibrium between the accre- 
tion and annihilation of WIMPs will be very roughly of 
the order of a year and the annihilation rate of WIMPs in 
the star will quickly reach a steady state such that equa- 
tion ([4]) will be equal to zero and all of the DM incident 
on the star can be assumed to go instantaneously into 
increasing the luminosity of the star. The luminosity of 
the star purely due to the accretion and subsequent anni- 
hilation of DM will therefore be given by (3t> 2 ,, c ^> 2(v 2 ) 
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which can be regarded as the minimum expected lumi- 
nosity of a star in an environment of weakly interacting, 
self annihilating DM of density pdm and velocity v. 

Observations of faint white dwarfs have been made in 
the central region of globular cluster M4 using the equiv- 
alent of UVI filters on the Wide Field/Planetary Camera 
(WFPC2) of the Hubble space telescope Observa- 
tions were taken in three fields, the closest one to the 
centre of the cluster with the planetary camera centered 
at 0.5r c ~ 0.2pc. 

We assume that the White dwarfs are perfect black 
body emitters, and we use their V-I magnitudes to cal- 
culate their temperature. We then use this temperature, 
plus the distance to globular cluster M4, to calculate the 
absolute bolometric magnitude of the stars, and conse- 
quently their luminosity. The resulting luminosities and 
temperatures of the white dwarfs are plotted in figure [2] 
In the same field, there will be stars further away from 
the centre of the cluster than 0.2 pc. However, due to 
the centrally peaked nature of the King profile, we calcu- 
late that more than 99% of the stars along a line of site 
at angular distance from the centre of the cluster corre- 
sponding to 0.2 pc will be at radii less than 5 pc from the 
centre of the cluster and will therefore be surrounded by 
the same density of DM. 

We derive the radii of the white dwarfs using their 
temperatures and luminosities and the assumption that 
they (and the sun) are black-body emitters, i.e. 
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which fits well with the radius expected for 0.5M Q white 
dwarfs using the standard mass-radius relations for white 
dwarfs. 

In figure [2] we have also plotted two lines correspond- 
ing to the luminosity one would expect due to the accre- 
tion and subsequent annihilation of WIMPs if the white 
dwarfs were placed in a region where the DM density is 
50 Mq pc -3 , as we expect in the inner region of the glob- 
ular cluster where they are located. We plot lines rather 
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FIG. 2: Luminosity vs. temperature for the white dwarfs in 
the inner field of the observations made in paper [l8| ] (data 
points). Also plotted (lines) is the minimum luminosity ex- 
pected for white dwarfs radiating as black bodies, for WIMP- 
nucleon cross sections of 10 -43 cm 2 and 10 -44 cm 2 . See text 
for details. 



than points because we calculate the expected luminos- 
ity and temperature for white dwarfs of different masses, 
and hence radii. The two lines correspond to two differ- 
ent WfMP-nucleon cross sections, IO -44 cm 2 and 10~ 43 
cm 2 , ff such a density of WIMP DM did exist in the 
centre of globular clusters with a WIMP-nucleon cross 
section corresponding to one or other of the lines, then 
we would not expect to see any white dwarfs lying below 
the line. 

Since the density of DM is robust to the total amount 
of matter in the cluster, this analysis suggests that if a 
core of DM remains in globular clusters then we are able 
to place a constraint upon the WIMP nucleon cross sec- 
tion between 10 -43 and 10 -44 cm 2 which is very compet- 
itive with the latest experiments. An alternative, more 
conservative statement would be that it is possible with 
this analysis to constrain the combination of WIMP nu- 
cleon cross section and DM density in globular clusters. 
It would therefore be interesting if WIMPs were detected 
in the laboratory with cross sections close to this level 
10~ 44 cm 2 , as it could open a new branch of DM astron- 
omy using white dwarfs. 

Since 1997, the same group who obtained the data used 
in this paper have made much more detailed studies of 
NGC 6397 but at a much larger radius from the centre 
of the cluster. Their observations are rather better than 
those used in this paper, enabling them to observe the 
super cool white dwarfs becoming more blue as H2 forms 



in their atmosphere, blocking the emission of redder fre- 
quencies. By obtaining colour-magnitude diagrams of 
this detail at different radii in the same cluster and com- 
paring them we predict that it should be possible to sig- 
nificantly increase the sensitivity of these constraints. 



IV. CAPTURE OF DM ONTO NEUTRON 
STARS 

In this section we will look at the accretion of DM par- 
ticles onto compact objects and their subsequent effects. 
As in the case of white dwarfs, the time scale for equilib- 
rium between capture and annihilation will be very short 
compared to typical astrophysical timescales. 

The most obvious possible observable signal of the ac- 
cretion of DM onto neutron stars is also the simplest 
to calculate, namely the heating of the star, with con- 
sequent increase of the surface temperature, due to the 
energy injected by annihilating DM particles. We follow 
the same procedure as that outlined in the previous sec- 
tion for the case of white dwarfs. For neutron stars and 
WIMP-nucleon cross sections close to the experimental 
limit, the surface area of the star is larger than the sum 
of the cross sections for the individual nuclei. (Unlike 
a nucleus with a large atomic number, a medium with 
a constant number density of nucleons gives rise no res- 
onant enhancement, although the presence of rod and 
sheet like structures in the neutron matter may give rise 
to an order of magnitude enhancement as they can for 
the neutrino- nucleon cross section [40|.) 

The expected heating of neutron stars is shown in fig- 
ure [3] for two different density profiles, namely r _1 and 
r -15 , both normalised so that the density at the solar 
radius 8.5 kpc from the centre of the galaxy is 0.3 GeV 
cm~ 3 . We have assumed a cross section large enough 
so that the surface area of the star is the effective cross 
section. We have also taken the most basic approach, 
namely that the neutron stars are in steady state and 
instantly emit the luminosity injected as black body ra- 
diation, (for a more detailed recent treatment, see fill ]) 
This figure shows that only at the solar radius can one 
expect the heating due to DM to be less than the min- 
imum temperature one might expect for a neutron star 
which started at some high temperature 13 billion years 
ago (see estimate in previous section.) It is however dif- 
ficult to observe the temperature of cool neutron stars 
- the small radius of these objects means that only ex- 
tremely hot neutron stars can be detected through their 
black body emission. It might be possible to place con- 
straints on the heating due to DM by looking at the non- 
observation of black body temperature from the surface 
of neutron stars which have been identified because they 
are pulsars. 

The CHANDRA x-ray telescope has also observed 
many hot compact objects in the central parsecs of the 
galaxy , although any cooler neutron stars in the same 
region would presumably be undetectable due to the dust 
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FIG. 3: Temperature of neutron stars due to the accretion of 
WIMPs as a function of distance from the centre of the galaxy. 
The two lines correspond to two different density profiles, p oc 



r and p oc r 



in the plane of the galaxy which effectively cuts out the 
transmission of photons below 2 keV down to energies less 
than the visible spectrum. The heating effect outlined 
here is relatively large though and it should be possible 
to find some constraining situations. 

Now we will consider other effects of the build up of 
DM in neutron stars. The reader should be aware of 
relevant previous work which considers the build up of 
charged DM in neutron stars [43| . 

The timescale for DM inside the star to thermalise with 
the background neutrons can be estimated - DM parti- 
cles falling onto a neutron star will be semi-relativistic 
and therefore need to lose an amount of kinetic energy 
roughly equivalent to their own mass through collisions 
with nuclei before they reach equilibrium. Kinematics 
show that the typical energy exchanged when the DM 



is more massive than nuclei will be AE ~ M„ 



I211 



so that as the DM moves through the neutron matter at 
speed close to c the timescale for it to lose energy is given 
by 
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where a is the DM nucleon cross section and n is the 
number density of nucleons. For example, if we would 
like a thermalisation time scale of less than a million 
years, we would require that the DM nucleon cross sec- 
tion a > 10 _60 (mdm/'Ti nuc )cm 2 . Once thermalised, the 
thermal radius within which the majority of thermalised 



DM particles will be located is given by 
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where we have assumed that the phase space density 
of DM particles is low enough that Maxwell-Boltzmann 
statistics are still valid, an issue we will return to below. 
The number of particles required to reach equilibrium 
between annihilation and capture is then given by 
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where a ann v = 10~ 26 cm 3 s~ 1 is the appropriate value for 
a thermal relic that achieve . The capture rate of 10 29 
s _1 corresponds to a rather large density of DM, but one 
which could be feasibly found at the galactic centre for 
100 GeV WIMPs. Throughout the rest of this paper, we 
will assume an accretion rate onto the neutron star of 
T c = 10 29 (100GeV/m dm )s- 1 . 

As the density of DM rises in the core of the neutron 
star, there will come a point at which the self-gravity of 
the DM core is greater than the gravity due to the bary- 
onic matter within the same volume. We label the num- 
ber of WIMPs in the star at the moment when this oc- 
curs with Ngo, defined by 47rr 3 /l p c /3 = NsG m dm- This 
means that if N eq > Nsg the core of DM will become 
self gravitating. For a given DM particle mass, this will 
occur when there are 
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DM particles in the star. The region of the parameter 
space where this happens in shown in fig. |4j 

This result only applies in cases where the Maxwell- 
Boltzmann statistics is relevant (assuming the DM can- 
didates are fermions). Consider a ball of radius R 
made up of N non-relativistic degenerate fermions of 
mass m and temperature T gravitating (not necessarily 
self-gravitating) due to a matter density p. These pa- 
rameters will be related by the approximate expression 
R 7 p 2 Gm - N 5/3 . The distribution of DM particles will 
form a degenerate 'dark star' when two conditions are 
fulfilled, the first is that the temperature is less than the 
Fermi energy, T < n 2 / 3 /mdm which occurs for a number 
of DM particles greater than 



Fermi 



= 3 x 10 



37 



T V 

ViqskJ 



gem 



Pc 



3/2 



(18) 
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independent of the DM parameters. The second require- 
ment is that the radius R is larger than the thermal ra- 
dius r t h, a situation which occurs when the number of 
DM particles N is larger than 



N d 



2.5 x 10 



\10 5 K J 



21/10 



/i0 14 gcm- 3 \ 9/1 ° /100GeV\ 3/2 



(19) 



Newly accreted DM particles would fall down onto the 
surface of this degenerate dark star, making it grow in 
mass and shrink in radius. For DM candidates above few 
hundred GeV, and adopting typical NS parameters, DM 
particles become self-gravitating before these degeneracy 
conditions are met. We also note that for a typical WIMP 
candidate with mass close to 100 GeV and the cosmolog- 
ically favoured DM self-annihilation cross section, we do 
not expect the DM to form such a degenerate dark star or 
a self-gravitating core. However, more massive particles, 
especially if strongly-interacting with ordinary matter, 
would be sev erely constrained, like for example SIMPzilla 
particles [H, |45| | which are a version of superheavy DM 
which could be produced non-thermally during^ inflation 
and may have masses as large as 10 12 GeV [4j, l5fl. These 
particles can interact with standard model matter with 
a cross section as large as a ~ 10 _24 cm 2 while their self 
annihilation cross section is limited by unitarity to be 
less than (a ann v) — msiMPzilla 0, Ea|- Note also that 
indirect detection arguments have been used recently to 
set an upper limit on the annihilation cross-section which 
is stronger than the unitarity bound for relatively light 
DM candidates 46] . The amount of material along a path 
through the centre of a neutron star will correspond to 
10 44 cm 2 so while each collision will only allow a maxi- 
mum momentum exchange of Q ~ m neutr on there will be 
enough material to slow down and stop these particles, 
so that the geometrical cross section will be appropriate 
in equation |T]) . The capture rate in regions of the galaxy 
with DM densities similar to that in the solar system will 
therefore be T c ~ 10 12 s~ 1 . Since the number of such high 
mass particles in a neutron star required to form a self 
gravitating core will be of the order of N$g = 10 16 , such 
a core would quickly form in a matter of hours. If the 
heavy DM candidate has a self annihilation cross section 
anywhere near the unitarity bound it will then instantly 
annihilate with itself. A similar fate would await other 
heavy DM candidates although, as discussed above, DM 
candidates which are too light will feel their own degen- 
eracy pressure before they collapse. 

We can estimate the maximum energy which could be 
released in such a burst of annihilations due to the self- 
gravitating collapse of a cloud of DM particles. Equa- 
tions (fT7| and HU) show us that the biggest release of 
energy in this way would be the sudden collapse and self 
annihilation of around 10 36 DM particles of mass 10 TeV. 
Such an injection of energy could in principle instigate a 
QCD phase transition in the central few centimetres of 



the neutron star, changing the equation of state in that 
region and possibly leading to a change in the observed 
rotation of pulsars. While glitches are observed in the pe- 
riodicity of pulsars [47] , more detailed study would have 
to be done to quantify the magnitude and observability 
of any such effect. 

The equation of state of any degenerate core will of 
course become relativistic when the number of particles 
reaches the Chandrasekhar limit, which corresponds to 
the number 



N C ha 



\m dm 



1Q51 ^lOOGeVy (20) 



V m dm J 



which for a 100 GeV DM particle is roughly to the mass 
of a planet like Mars. For a SIMPzilla candidate with 
the properties outlined above, the Chandrasekhar limit is 
much less, of the order of 1000 tonnes. If that is the case, 
the accretion of matter onto the newly formed mini-BH 
would rapidly destroy the star over very short timescales 
(see Ref. [43j] for a detailed discussion of the growth of the 
mini black- hole, and consequent destruction of the NS). 

We have summarised the different behaviours in figure 
HI where we have assumed a large DM accretion rate of 
T c = (100GeV/m dm )10 29 s- 1 . If the DM densities advo- 
cated in this paper can be reached in nature, the obser- 
vation of compact objects in those regions would place 
strong constraints on the particle physics parameters of 
DM particles. 



V. CONCLUSIONS 

In this work, we have investigated the observational 
consequences of the capture of DM onto astrophysical 
compact objects. Although white dwarfs and neutron 
stars are relatively hot, making it a challenge to detect 
any heating of such objects due to DM annihilations in 
their cores, we have presented one situation, namely the 
possible accretion of WIMPs onto cool white dwarf stars 
in globular clusters, where we argue it is feasible that such 
heating could be detected. The presence or otherwise of 
DM in the core of globular clusters is model dependent, 
but we have argued that if DM does exist in the cluster, 
then the density at the location of the white dwarfs in 
question will not be a strong function of the total amount. 

We have also estimated the heating of neutron stars 
due to the accretion of WIMP DM (see also the re- 
cent [4lj). Finally, we have explored the parameter space 
of DM mass and self-annihilation cross section, focusing 
on the accretion of such particles onto neutron stars in a 
region of the galaxy where the density of DM is high. We 
have found that there are a number of different outcomes 
of DM accretion, depending upon the exact parameters 
chosen, including the creation of degenerate dark stars in- 
side the neutron star, or self-gravitating cores. It would 
be interesting to study the fate of such configurations 
when they become unstable and collapse, in order to look 
for observable consequences. 
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FIG. 4: Different outcomes of the accumulation of DM inside a neutron star, in the DM mass vs. annihilation cross section plane, 
for a capture rate of T c = (100GeV/m<im)10 29 s _1 . In the top left corner we show the results of a scan of the supersymmetric 
parameter space, where neutralino models are compatible with accelerator and cosmological constraints (as obtained with 
DarkSUSY 48]). The solid line is relative to viable DM models in Universal Extra-dimensions (see text for further details). 
In the shaded region below the dashed line, DM particles become self-gravitating before equilibrium between capture and 
annihilation is reached. For models on the right of the dotted vertical line, particles reach the critical mass for gravitational 
collapse in less than 1 Gyr. These models can likely be ruled out, since they lead either to large injection of energy in the NS 
core, or to gravitational collapse to a Black Hole, rapidly destroying the star. 
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